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Abstract. In this paper, we wish to point out the possibility of using a complex 
scalar field to account for the inflationary stage and the current acceleration. 
By the analysis of the dynamical system and numerical work, we show that the 
amplitude of the complex scalar field plays the role of the inflaton whereas the 
phase is the quintessence field. The numerical solutions describe heteroclinic 
orbits, which interpolate between an unstable critical point and a late-time de 
Sitter attractor. Therefore, this model is more natural to explain the two stages 
of acceleration. 
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1. Introduction 

Astronomical observation on the cosmic microwave background(CMB) anisotropy 
supernova type la(SNIa) and SLOAN Digital Sky Survey(SDSS) converge on 
that our Universe is spatially flat, with about 70% of the total density resulting from 
dark energy that has an equation of state w < —1/3 and drives the accelerating 
expansion of the Universe which began at a redshift of order one-half. The origin 
of the dark energy remains elusive from the point of view of general relativity and 
standard particle physics. To represent dark energy, several candidates have been 
suggested and confronted with observation: a cosmological constant, quintessence 
with a single field 0] or with N coupled field a phantom field with a canonical 
|B] or Born-Infield type Lagrangian 0, fc-essence |S] and the generalized Chaplygin 
gas(GCG) [n]. Among these models, the most typical ones are cosmological constant 
and quintessence which has caught much attention ever since its invention. 

The idea of inflation is legitimately regarded as an great advancement of modern 
cosmology: it solves the horizon, flatness and monopole problem, and it provides a 
mechanism for the generation of density perturbations needed to seed the formation 
of structures in the universe ^Uj- In standard inflationary models the physics 
lies in the inflation potential. The underlying dynamics is simply that of a single 
scalar field rolling in its potential. This scenario is generically referred to as chaotic 
inflation in reference to its choice of initial conditions. This picture is widely favored 
because of its simplicity and has received by far the most attention to date. The 
properties of inflationary models are also tightly constraint by the recent result from 
the observation. The standard inflationary ACDM model provides a good fit to the 
observed cosmic microwave background (CMB) anisotropy. Peebles and Vilenkin 
proposed and quantitatively analyzed the intriguing idea |12j that a substantial 
fraction of the present cosmic energy density could reside in the vacuum potential 
energy of the scalar field responsible for inflation (quintessential inflation). After 
that, there were some models to be presented in succession. 

In this paper, we wish to point out the possibility of using a complex scalar 
field to account for the inflationary stage and the current acceleration. We put the 
emphasis on the study of the dynamics. According to the phase space analysis and 
numerical calculation, we show that how the amplitude of complex field plays the 
role of the inflaton whereas the phase is the quintessence field. There are stable and 
unstable critical points for our model. In our case, one of the unstable critical point 
is corresponding to the inflationary stage at very first moments of the Universe and 
the stable critical point is corresponding to the second stage of accelerated expansion 
began at a redshift z ~ 1.5 and is still operative, which alleviates the fine tuning 
problem. The equation of state w varies with the cosmic evolution and approaches 
towards -1 asymptotically showing the existence of a cosmological constant at late 
times. The heteroclinic orbits 15 connect unstable and stable critical points. The 
numerical calculation shows that brought the Universe back to the usual Friedman- 
Robertson- Walker expansion, then the second stage of accelerated expansion began at 
z ~ 1.5. 
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Since current observations favour a flat Universe, we will work in the spatially flat 
Robertson- Walker metric, 



ds^ = dt^ - a^{t)d-x^, 
the Lagrangian density for the spatially homogeneous complex scalar field $ is 

= {g'" (5^$*) (5.$) - vm) - c/(0), 

when consider the presence of barotropic fluid, the action for the model is 

1 
2^2 



S = / d'^x^ 



(1) 



(2) 



(3) 



where g is the determinant of the metric tensor g^i,, k} = SttG, TZ is the Ricci scalar, 
and p-y is the density of the fluid with a barotropic equation of state — (7 — l)^^, 
where < 7 < 2 is a constant that relates to the equation of state by w = 7 — 1. 
Writing the complex field $ as (see 1131) 



(4) 



Substituting |0J into Q and varying the action, one can obtain the Einstein equations 
and the equations of motion for the scalar field as 
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1 



a2 



= -3H{py+p^), 



P 



'+\3H + 2^] e+'—u'[e) = 0, 



(5) 

(6) 
(7) 

(8) 

(9) 



where H is the Hubble parameter, dot and prime denote derivatives with respect to t 



and (j) (or 9) respectively. We define = \4P' + 1/(0) and pg = 57^^^ + U{e). 



3. Phase space 

In this section, we investigate the global structure of the dynamical system via 
phase plane analysis and compute the cosmological evolution by numerical analysis. 
Introduce the following dimensionless variables: 



= -7^<p, yi - 



V6H 



V3H ' 



^2 = -—p^O, ya - 



f Vgh 



V3H ' 
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^ /tJ7(6')' ^ C/'2(6i) ' 
/r 

0=—, iV = lna, (10) 
and dimensionless constant 

^ = fK. (11) 
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(12) 



Note that / = 0{Mpi), n ~ \/87rG, G = m~j^f (natural units) and mpi — \/8ttMpi, 
so we could have ^ = 0{1). It is useful when we do numerical calculation. 
Now, the equations ©-© become the following system: 

dxi 
dN 

dX2 

dN 
dyi 
dN 
dy2 
dN 
dz 
dN 

^ = -^^A?Xl(^l-l), 

^=-XUzX2{T2-l). 

Also, we have a constraint equation 

+ + ^ = (13) 

where 

n,^!f^^xl+yl (14) 

ne^^^xl + yl (15) 

The equation of state for the complex scalar field could be expressed in terms of the 
new variables as 

^ ^ xj + xl-yl- yl 
xj+xl+yf + yl 

According to the definitions, the parameters Ai, Fi, A2, r2 become constants and 
are equal to Ai(o), 1, A2(o) and 1, respectively. Under the circumstance, the equations 
p2|l constitute an autonomous system as follows, 

'^■^i 3 r 2 2 2 2\ , r./ 2 , 2\n r, , 2 , / 3 X 2 



dN 2^^ ^^^^ -x\~xl-yl- yl) + 2{x\ + xD] - 3xi + xlz + i^^mvl, 
^ = ^X2 [7(1 -xl-xl-yl- yl) + 2{xl + xj)] ~ 3x2 - X1X2Z + ^X2{o)^ylz, 



A cosmological model with complex scalar field 



5 



^ = ^2/1 [7(1 - a;? - - 2/1 - vl) + 2(a:? + 2^2)] - ^^^m^iyi, (17) 

^ = [7(1 -xj-xl-yf- yl) + 2(a;^ + xD] - ^A2(o)C2;2y22, 

dz 2 

= —x-\ z . 

dN 

When barotropic matter is under consideration or the equations of motion are 
too difficult to solve analytically, phase space methods become particularly useful, 
because numerical solutions with random initial conditions usually do not expose all 
the interesting properties. In table 1, we list the critical points and the cosmological 
parameters there. To gain some insight into the property of the critical points, 
similarly as in ^1], we investigate the stability of these critical points. For the critical 
points listed in table 1, we find the eigenvalues of the linear perturbation matrix, see 
table 2. For stability we require the all 5 eigenvalues to be negative. 

The corresponding conclusions of numerical calculation are shown in Fig^ and 
FigEl This numerical solutions describe heteroclinic orbits, which interpolate between 
an unstable critical point (case i) and a late-time de Sitter attractor (case iii). 
Therefore, this model is more natural to explain the two stages of acceleration. A 
point worth emphasizing is that the ordinary matter (radiation and dust) affect the 
evolution of the scalar field via their contribution to the general expansion of the 
Universe because of the couples can be neglected in the model. In FigQ] the behavior 
of the equation of state parameter is shown. At initial time, u; ~ — 1 (we choose 
A 1(0) = 0-5, then w ~ 0.916) corresponds the inflationary phase, then it increases and 
becomes positive. After arriving at the value 1/3, the Universe comes to the radiation 
dominated epoch and w stays on a broad platform. Next, w drops to zero and stays 
on a narrow platform. Finally, w drops below zero and approaches to —1, which 
corresponds second stage of accelerated expansion. The evolution of cosmic density 
parameters are shown in Fig[21 The phase part contribution VIq stays at i7 = at the 
very first moments of the Universe and becomes 1 in late-time, which plays the role 
of quintessence field. On the other hand, the amplitude part contribution fi^ plays 
the role of inflaton. During the whole evolution, the radiation energy density and 
the dust matter energy density become dominate respectively. Therefore, we see that 
the constraints arising from cosmological nucleosynthesis and structure formation are 
satisfied. 



Case Critical points {xi, X2, yi, 2/2, z) ^ct> w 



(i) 


v/%,0, v'i-%,0,0 


1 





_1 + ^ 

^ 3 


(ii) 




^l(O) 
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7-1 


(iii) 


0, 0, 1, 0, 





1 


-1 


(iv) 


±1, 0, 0, 0, 


1 





1 


(v) 


0, ±1, 0, 0, 





1 


1 


(vi) 


0, 0, 0, 0, any 









Table 1. The properties of the critical points. 
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case corresponding eigenvalues stability 

~(i) 0, -3 + -3 + %i, -37 + A2(o) unstable 

(ii) 0, f, 1(7-2), 

3{A^(0)(7-2)-y^A^(0,(7-2)[A^^,,^(97-2)-247'^]} 



'^l(O) 

(iii) -3, 0, 0, 0, -37 stable 

(iv) 3, 0, 0, 6 - 37, 3 T ^1^1(0) unstable 

(v) 3, 3, 0, 0, 6 - 87 unstable 

(vi) 0, ^, ^, -3+^, -3+'^ unstable 



Table 2. The eigenvalues of the critical points. 
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Figure 1. Evolution of the equation of state at the presence of radiation and 
matter. The curves correspond to the different initial conditions. 
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Figure 2. Evolution of cosmic density parameter of quintessence energy density 
Q,e, inflaton energy density radiation energy density Hr and matter energy 
density Clm- The curves correspond to the different initial conditions. 
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4. Conclusions 

In this paper, we liave discussed the cosmological implication of a complex scalar model 
of quintessence inflation with a barotropic fluid. Using the numerical calculation for 
this model, we show that the phase contribution is negligible at the early epoch of the 
universe while it becomes dominate with the time evolving and the evolution of the 
amplitude part is the opposite, which is a viable way to unify the description of the 
inflationary stage and the current accelerated expansion. Analysis to the dynamical 
evolution of the complex scalar model indicates that it admits a late time attractor 
solution, at which the field behaves as a cosmological constant. Obviously, attractor 
and heteroclinic orbit are both alleviate the fine tuning problem. 
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